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CHAPTER  I 


THE  PROBLEM 


Introduction 


In  recent  months,  the  United  States  has  decided  to  deploy  an  anti -ballistic 
missile  system.    This  system  employs  two  missiles,  one  for  high  altitude,  long 
range  interception  (the  Spartan)  and  the  other  for  low  altitude,  short  range  inter- 
ception (the  Sprint) .    The  mission  of  the  Sprint  is  to  provide  "insurance",  i.e. ,  to 

get  the  targets  that  Spartan  missed.    Sprint  will  intercept  targets  at  altitudes  of 

*     1 
100,000    feet  or  less,  and  will  have  an  acceleration  capability  of  over  lOOg  (11) . 

The  difficulty  of  the  missile  interception  problem  is  primarily  due  to  the  short 

time  available  for  interception.    An  ICBM  warhead  just  entering  the  atmosphere  is 

traveling  at  nearly  18,000  miles  per  hour.    Thus,  every  second  of  fli^t  brings  it 

five  miles  closer  to  the  target.    The  defending  radar  system  has  very  little  time  in 

which  to  detect  the  warhead  and  define  its  trajectory  —  less  than  three  minutes, 

according  to  one  source  (10) .    It  is  advantageous  to  fire  the  interceptor  as  late  as 

possible,  in  order  to  give  the  radar  maximum  time  to  track  the  target  and  thus  to 

refine  the  trajectory  calculations .    The  only  optimum  trajectory  for  the  interceptor 

is  a  minimum -time  trajectory,  because  this  will  allow  the  interceptor  to  be  fired 

* 

as  late  as  possible  (9) . 

■  ,       '  / 

I 

/ 
/ 

^Numbers  in  starred  parentheses  stand  for  the  numbered  references  at  the  end 
of  this  report. 


From  another  point  of  view,  it  should  be  noted  that,  for  a  constant-thrust  in- 
terceptor, minimum  time  imlplies  minimum  fuel  consumption,  so  the  interception 
could  be  made  at  the  largest  possible  distance  from  the  launch  point. 

With  the  foregoing  facts  in  mind,  the  purpose  of  this  report  was  twofold: 

1)  To  derive  the  necessary  conditions  for  a  minimum -time  trajectory  to  a 
fixed  point  in  range -altitude  space, 

and 

2)  to  numerically  integrate  the  equations  derived  in  part  one  for  a  hypo- 
thetical "Sprint -class"  interceptor  and  use  the  results  in  a  realistic 
scheme  for  interception.  ^ 

physical  Assumptions 
A  picture  of  the  missile  in  flight,  along  with  the  forces  acting  on  it,  is  pre- 
sented in  Figure  1,  and  Table  1  is  a  list  of  nomenclature.    The  assumptions  made 
in  this  development  were  : 

1)  Planar  motion 

2)  No  lift ;  thus  the  only  way  to  control  the  fli^t  path  is  by  thrust  vectoring 

3)  Constant  thrust 

4)  Flat  earth,  constant  g 

5)  Drag  is  a  function  of  altitude  and  velocity,  i.e. ,  D  ■   D(h,V) 

1 

6)  Mass  is  a  linear  function  of  time,  i.e. ,  m  ■  hIq  -  m't. 

The  problem  is  to  find  the  time  histories  of  the  problem  variables  x,  h,  V,  A,  and 
B  such  that  the  time  cf  flight  between  (xi.hj)  and  (xf.hf)  is  a  minimum.    . 


-J  I 

Prime  superscripts  denote  differentiation  with  respect  to  time.  ' 


TABLE  I 


Nomenclature 


T Thrust  magnitude,  Ibf. 

D i Drag  magnitude,  Ibf. 

2 
g Acceleration  of  gravity,  32. 174  ft/sec 

m Mass  of  Interceptor,  slugs 

f  subscript Final  value 

1   subscript •  •  •    Initial  value 

"Problem  Variables'.'  ' 

X Range,  feet 

h •  Altitude,  feet 

V Velocity  magnitude,  ft. /sec. 

A  .  • •  •  Angle  between  velocity  vector  and 

horizontal,  measured  from  horizontal, 
positive  counterclockwise. 

B .    Angle  between  thrust  vector  and  velocity 

vector,  measured  from  velocity  vector, 
positive  counterclockwise  ;  also  called 
"'control  variable"  and  "thrust  angle" 

"L  Variables";  .      , 

Lj    ......................  Lagrange  Multiplier  associated  with  x 

L2 "  "  "  "    h 

L3   .  . .        "  ••  "  "V 

L4 "  '•'  "  "A 


CHAPTER  II  i 


NECESSARY  CONDITIONS  FOR  A  MINIMUM-TIME  TRAJECTORY 

The  problem  of  finding  a  minimum  -time  trajectory  between  two  points  in 
range -altitude  space  resembles  the  classical  Brachistochrone  problem  of  the  Cal- 
cuius  of  Variations    ( 1 )    ,  and  a  similar  indirect  method  was  used  in  this  report. 
Specifically,  the  problem  was  formulated  as  a  Mayer  problem  in  the  Calculus  of 
Variations.    A  short  outline  of  the  Mayer  problem  as  a  special  case  of  the  Bolza 
problem  is  given  in  Appendix  1. 

Equations  of  Motion 

From  Figure  1,  a  force  balance  on  the  missile  provides  the  following  equa- 
tions: ... 
mV  =    TcosB  -  mgsinA  -  D  (1) 
mVA'=    TsinB  -  mgcosA      •  (2) 
The  kinematic  relationships  are  -                                    * 

x'»    VcosA  '                            .                                  (3) 

h'=    VsinA  (4) 

Formulation  of  the  Variational  Problem 

There  are  many  ways  to  formulate  a  problem  to  minimize  time.    First  of  all, 
an  expression  for  the  time  of  flight  must  be  found.    Any  one  of  equations  (1)  to  (4) 


could  be  solved  for  dt  and  an  expression  for  the  time  of  flight,  formed.    From  (1) , 

tf  -  tj  -    l(m/(TcosB  -  mgsinA  -  D))dV  ; 

I 
from  (2) ,  I 

M  •  ■  ■     : 

tf  -  tj  -   j  (mV/(TsinB  -  mgcosA))dA  ;  | 


j  (ti 


from  (3) , 


tf  -  tj  -   I  (1/VcosA)  dx  ; 


from  (4) 

tf  -  ti"  I  (l/VsinA)dh. 
Any  of  these  formulations  gives  rise  to  a  Lagrange  problem,  that  is,  the  extremiza- 


tlon  of  an  integral  subject  to  side  constraints.    However,  since  the  problem  at  hand 
is  in  range -altitude  space,  it  would  seem  wise  to  choose  an, expression  with  x  or  h 
as  the  independent  variable  as  in  the  last  two.    At  least  one  author  has  formulated 
the  problem  in  this  manner  (7) . 

All  of  the  foregoing  methods  effectively  eliminate  time  from  the  calculations,    . 
leaving  velocity,  range    altitude,  or  path  angle  as  the  independent  variable.     A 
much  simpler  formulation  is  that  of  Mayer,  in  which  case  time  is  the  independent 

variable  ;  by  this  method  of  formulation,  a  more  general  class  of  problems  can  be 

* 

treated  with  the  same  equations  (3) . 


Mayer  Problem 


The  Mayer  problem  is  a  special  case  of  the  problem  of  Bolza,  and  both  are 
outlined  in  Appendix  1.    Applying  the  Mayer  formulation  to  the  present  problem, 
■we  must  first  rearrange  Equations  (1)  through  (4)  and  introduce  the  new  functions. 

Pi:  ''y 

Pl«»  x'  -  VcosA  "0,  '    (5)      j: 

P2»    h'  -  VsinA  =  0,  (6)      I 

,   ■■  .  i 

P3=    mV  -  TcosB  +  mgsinA  +  D  -  0,  (7)     '[.-. 

■■ . .  '  -  (. 

and  i.  ' 


P4-    mVA'  -  TsinB+mgcosA=    0. 
The  augmented  function,  F,  is  defined  as 

F-   L^Pi  +  L2P2  +  L3P3  +  L4P4 
where  the  L's  are,  in  general,  functions  of  t. 

feuler- Lagrange  Equations 

The  Euler- Lagrange  equations  for  this  problem  are  : 

'     l'i-o, 

t 

L3Dh  -  L'2  -0, 

L1COSA+  L2sinA  -  L3(Dv  -  m')  -  L4mA'  +  L*3m  "  0, 

LiVsinA  -  L2VCOSA  +  L3mgcosA  -  L4(mgsinA  +  mV  +inV) 

and 


0, 


(8) 


(9) 


(10) 
(11) 
(12) 
(13) 


LsTsinB  -  L4TCOSB  »  0. 


(14> 


Since  Equation  (14)  is  algebraic,  it  enables  us  to  eliminate  one  of  the  Vs,  spec- 

iiScally  L4  : 

L4  «   L3tanB.  (15) 

Upon  substituting  (15)  into  (12)  and  (13) ,  rearranging  the  equations,  including 

Equations  (5)  through  (8) ,  and  noting  that  Equation  (10)  merely  says  that  Li» 

constant,  we  arrive  at  a  set  of  seven  first-order  ordinary  differential  equations 

in  the  seven  unknowns  x,h,V,A,B,L2,  and  L3  :        . 

'•        x'-    VcosA,  .  (16) 

h'-    VsinA,  '  (17) 

V-     l/m(TcosB  -  mgsinA  -  D) ,  (18) 

A'»     l/mV(TsinB  -  mgcosA) ,  ^    ■  (19) 

L^-   LsDh.  (20) 

L'3-    l/ni(L3(I\  +  mA'tanB  -  m')  -  LjcosA  -  L2SinA) ,  (21) 

i  ■ 
and 

B'  -    (cos^B/mVL3)(LiVsinA  -  L2VCOSA+  L3(mgcosA  -  mgsinAtanB- 

mV'tanB  -  m'VtanB)  -  mVL'3tanB).  .  (22) 

This  set  of  equations  appears  to  be  of  the  seventh  order  in  the  aggregate,  and  one 

would  normally  expect  to  need  seven  boxindary  conditions.    The  question  which 

then  arises  is  "how  many  conditions  can  be  imposed  on  the  problem  variables 

(essential  boundary  conditions)  and  how  many  are  natural  end  conditions  involving 

the  L's?"    In  short,  how  many  physical  end  conditions  can  be  specified  without 

over -determining  the  problem?   The  answer  is  provided  by  an  analysis  due  to 

KeUey  (3). 
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"Effective  Order"  of  a  Set  of  Equations 

According  to  Kelley ,  in  a  Mayer  problem ,  the  number  of  boundary  conditions 
which  can  be  imposed  on  the  problem  variables  is  equal  to  the  "effective  order"  of 
the  set  of  equations ,  and  this  effective  order  is  always  one  less  than  the  apparent 
order.    This  is  due  to  the  linearity  and  homogeneity  of  the  equations  in  the  L  vari- 
ables, and  can  be  seen  in  the  following  way.    If  the  initial  (or  final)  values  of  all 
three  L's  are  multiplied  by  the  same  factor,  there  is  no  change  in  the  problem 
variables.    Only  the  ratios  of  the  L's  are  significant. 

Thus,  in  the  present  problem,  the  effective  order  of  the  set  of  equations  is 
six,  so  six  essential  boundary  conditions  can  be  imposed  on  the  problem  variables. 

Boundary  Conditions 

In  this  case,  the  essential  boundary  conditions  chosen  were  the  following  : 
Attj=    0:    x(0)-  0,    h(0)=    0,    V(0)-     Vj,  and  A(0)  »    Aj.  (23) 

At  tf:  X  (f)=     xf  and  h  (f)  «  hf .  (24) 

The  other  conditions  are  natural  end  conditions,  and  they  are  obtained  from  the 

transversality  condition. 

The  transversality  condition  for  the  present  problem  is  • 

((1  -  Lj  VcosA  -  L2VsinA  -  L3(TcosB  -  mgsinA  -  D)  -  LatanBCTsinB  -  I 

f  '.■'"•; 

mgcosA) )  dt  +  Lidx  +  L2dh  +  LamdV  +  LamVtanBdA)^-    0.  /     j. 

From  Equation  (23) ,  ! 

i    ■ 
dxj-    dhj-    dVi-    dAj-    0.  ]' 


From  Equation  (24) ,         ■ 

dxf  =  dhf  "0. 
Since  t^  is  fixed  at  zero,  , 

dq-    0.  ■  •'  ,■;.•■ 

Since  Vf  is  open,  dVf  is  arbitrary  ;  thus, 

L3f  mf    0, 
and  since  mf  /  0, 

L3f-    0.  (25) 

Note  that  Af  is  open,  but  Lgf  "  0  ;  therefore,  no  additional  information  is  gained 
from  this  fact. 

Since  time  is  to  be  minimized,  tf  cannot  be  fixed  ;  therefore,  dtf  is  arbitrary. 
Combining  this  fact  with  Equation  (25) ,  we  have 

(1  -  LjVcosA  -  L2VsinA)f  »  0, 
or,  using  Equations  (3)  and  (4), 

(Lix' +  L2h')f  -   1.  (26) 

This  exhausts  all  possibilities.    A  total  of  eig^t  boundary  conditions  have  thus  been 
obtained.  -  . 

Since  Ljis  constant  and  only  ratios  of  the  L's  are  significant,  Lj  can  be  set 

equal  to  1,  and  the  ratios  L2/LI  and  L3/L1  become  merely  L2  and  L3  respectively. 

Also,  condition  (26)  becomes  (x'  +  L2h')f  •■  1. 

Chapter  III  of  this  report  contains  a  discussion  of  the  application  of  these 

/'      .      . 
boundary  conditions. 


CHAPTER  III 
NUMERICAL  INTEGRATION  OF  EULER -LAGRANGE  EQUATIONS 

The  problem  as  defined  in  Chapter  11  involves  four  initial  conditions  and  four 
final  conditions.    If  an  analytical  solution  to  the  set  of  differential  equations  could 
be  found,  this  would  pose  no  problem.    However,  since  no  solution  could  be  found, 
it  was  necessary  to  use  a  numerical  integration  process,  integrating  from  a  spec- 
ified set  of  initial  conditions  and  trying  to  converge  on  the  desired  final  conditions 
by  varying  those  initial  conditions  which  were  not  essential  boundary  conditions. 

In  the  present  problem,  the  desired  final  conditions  are  specified  by  Equations 
(24) ,  (25) ,  and  (26) .    The  non-essential  initial  conditions  which  can  be  varied  to 
obtain  the  desired  final  conditions  are  Bq,  L20  andLso-    It  is  conceivable  that,  with 
enough  diligence  and  experience,  one  could  converge  on  the  four  specified  final 
values  at  once,  but  this  possibility  was  considered  to  be  nil  under  present  circuni- 

stances.    At  this  point,  then,  it  was  thought  advisable  to  employ  a  device  used  by 

» 
Vincent  (8)  and  called  a  "flooding  technique". 

Flooding  Technique 

This  technique  is  essentially  the  following  t  .  .' 

1)  The  essential  boundary  conditions  at  the  initial  point  are  imposed  ; 

2)  arbitrary  initial  values  of  B  ,    L2  and  L3  are  imposed  ; 

/ 

3)  the  equations  are  integrated  to  missile  burnout ; 

4)  a  check  is  made  to  see  if  the  natural  end  conditions  are  met  (no  check  is 


10 


u 


made  on  Xf  and  hf  ) ; 
5)     the  initial  values  of  B,  L2  and  L3  are  varied  in  a  systematic  fashion,  and 
the  process  is  repeated  until  the  natural  end  conditions  are  met;  Xf  and  hf 
are  allowed  to  fall  where  they  may  ;  in  this  way,  optimal  trajectories  to  a 
field  of  endpoints  may  be  found. 
It  is  important  to  note  that  this  process  is  carried  out  for  a  given  set  of  es- 
sential initial  conditions  —  that  is,  for  given  values  of  xj,  hf,  Vi  and  A^.     The 
endpoints  of  the  optimal  trajectories  found  by  this  process  form  a  "locus  of  optimal 
endpoints".    By  changing  one  or  more  of  the  essential  initial  conditions  and  perform- 
ing the  flooding  process  again,  a  different  locus  of  optimal  endpoints  may  be  found. 
In  the  present  investigation,  the  initial  values  of  x,  h,  and  V  were  left  at  a  fixed 
value  throughout ;    only  Aj  was  varied  over  a  range  of  values,  a  locus  of  optimal 
endpoints  being  obtained  for  each  value  of  Ai. 

Although  these  trajectories  are  planar,  the  concept  is  easily  extended  to  three- 
dimensional  space  by  simply  rotating  each  locus  around  the  H-axis.  A  "surface  of 
optimal  endpoints"  is  thus  generated  which  forms  an  "umbrella"  over  the  launch 
site.    As  noted  in  Chapter  II,  this  umbrella  represents  the  locus  of  points  farthest 
from  the  launch  site  at  which  interception  can  occur, 

'  Model  Interceptbr 

Before  numerical  results  can  be  obtained,  some  data  on  the  interceptor  must 

* 
be  given.    The  information  released  so  far  (10, 11)  regarding  the  Sprint  missile 

gives  only  diameter  (4, 5  feet) ,  length  (27  feet) ,  acceleration  capability  (over  100  g) , 
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and  interception  ceiling  (under  100,000  feet).    Therefore,  a  typical  "Sprint-type" 
missile  fitting  at  least  these  meager  specifications  was  proposed  in  this  report. 
Its  characteristics  are  given  in  Table  2.    The  missile  was  assumed  to  be  fired, 
like  the  Sprint,  from  a  launcher,  by  an  explosive  process,  with  a  fixed  initial 
velocity  of  500  feet  per  second.    This  launcher  is  assumed  to  be  of  such  design 
that  it  allows  the  missile  to  be  launched  with  an  initial  fligjat  path  angle  (Aj)  of 
anywhere  from  30  degrees  to  90  degrees. 

Information  regarding  the  atmospheric  density  profile  and  drag  function  used 
in  this  report  is  given  in  Appendix  2.    The  Runge-Kutta-Gill  integration  scheme  is 

discussed  in  Appendix  3. 

'.     ■  ■  ■  '        ,        ■     " 

I^Jumerical  Results, 

As  was  noted  previously,  the  initial  conditions  Bj,  L2i  and  L3i  were  to  be 
varied  in  an  attempt  to  converge  upon  the  final  values  L3f  «  0  and  j6  +  L«,h' ■    1. 
The  left-hand  side  of  the  latter  equation  was  set  equal  to  a  quantity  called  PARAM, 
and  the  goal  was  then  to  obtain  the  conditions  PARAMf  »  1  and  Lof  •»  0  simultan- 
eously. 

In  order  to  get  some  idea  of  the  influence  of  changes  in  Bj,  L2i  and  Lq-  on 
PARAMf  and  L3J,  parametric  studies  were  conducted(  2  ).    In  each  study,  the  initial 
value  of  one  of  the  variables  was  changed  while  the  initial  value  of  the  others  re- 
mained fixed .    This  was  done  for  values  of  Aj  of  30,  60,  and  90  degrees.     The 
results  of  these  studies  are  shown  in  Figures  2  to  6. 

From  these  graphs,  a  number  of  conclusions  were  drawn.    From  Figure  6 


TABLE  2 

SPECIFICATIONS  OF  MODEL  INTERCEPTOR 

■"'■•■  .  -       . 

Number  of  stages 1 

Launch  mass '. 850  slugs  (27,400  Ibm.) 

Burnout  mass .    100  slugs  (3,220  Ibm.) 

Propellant  Burning  rate 50  slugs/sec. 

Exhaust  velocity 10,000  fps. 

Thrust 500,000  Ibf.  (constant) 

Burning  time     15  seconds 

Thrust  acceleration From  26  to  156  g       '  . 

Frontal  area 20  square  feet  ^ 


/ 
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it  was  seen  that  B  has  very  little  influence  on  PARAM,  and  L    .    It  was  found,  how- 

1  I  oX 

ever,  that  B  strongly  Influenced  the  values  of  x,  and  h  .    Thus,  the  following  pro- 
cedure was  adopted  : 

1)  Set  B  "    0  and  converge  on  the  natural  end  conditions  by  varying  L 

i  2i       . 

and  L     only  ; 
3i 

2)  once  the  final  conditions  have  been  met,  vary  B.  to  obtain  a  "fan"  of 
near-optimal  trajectories  (in  the  sense  that  PARAM  and  L     are  almost 
at  their  required  values  )  to  a  field  of  x-h  endpoints  ; 

3)  again  vary  L     and  L     very  sligjitly  to  correct  the  slight  errors  in 


PARAMj  and  L„,  inducted  by  the  changing  B  . 


This  was  done  for  A   equal  to  30,40,50,60,70,80,  and  90  degrees,  and  the  resulting 

1 

"fans"  of  trajectories  are  shown  in  Figures  7  to  13.    B.  was  varied  over  a  range  of 
only  ±  10  degrees,  because  it  was  felt  that  anything  over  these  bounds  was  unrealis- 
tic from  a  physical  standpoint. 

For  each  trajectory  in  Figures  7  to  13,  the  desired  natural  end  conditions  were 
satisfied  as  nearly  as  was  possible  with  the  limited  time  and  experience  at  hand.    In 
all  cases,  the  value  of  PARAM   fell  within  the  limits  of  410,    This  was  found  to  be 
adequate  due  to  the  non-sensitivity  of  the  trajectory  parameters  to  small  changes  in 
PARAM..    However,  the  trajectory  parameters  were  sensitive  to  small  changes  in 
Lg,,  and  difficulty  was  experienced  in  approaching  the  desired  value  of  L~,  "0.    In 
all  cases,  the  value  of  Lo*  was  below  0.09,  and  in  most  cases,  it  was  below  0.05. 
Further  efforts  to  reduce  L_,  were  unsuccessful, 

oX 

In  the  remainder  of  this  report,  the  trajectories  in  Figures  7  to  13  are  called 
"optimal  trajectories"  ;  however,  they  are  only  near-optimal  at  best,  because  the 
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natural  end  conditions  were  not  precisely  satisfied. 

It  was  thought  that  one  reason  for  this  difficulty  was  that  no  bounds  were  placed 
on  the  control  variable,  B.    For  a  more  valid  physical  model,  limits  should  have 

I 

been  imposed  upon  B,  and  this  would  result  in  a  different  problem  ;  this  different 
problem  is  discussed  in  Appendix  4. 

Loci  of  Optimal  Endpoints 


In  the  present  problem,  the  missile  has  a  fixed  burning  time  of  15  seconds,  and 
all  of  the  optimal  trajectories  have  endpoints  with  a  final  time  of  15  seconds.    These 
endpoints  form  a  locus  of  optimal  endpoints,   or  "15  second  locus".    There  is  a  locus 
for  each  value  of  A.,  and  these  loci  are  shown  as  the  dotted  curves  at  the  end  of  each 
optimal  fan  in  Figures  7  to  13.    Only  five  trajectories  were  run  within  each  fan,  and 
the  loci  of  optimal  endpoints  were  found  mearedy  iy  drawing  a  smooth  curve  throu^ 
the  endpoints  of  the  five  trajectories  in  each  fan. 

In  Figure  14,  all  of  the  fans  of  optimal  trajectories  shown  in  Figures  7  to  13 
are  superimposed  on  one  graph.    It  can  be  seen  that  the  loci  overlap,  and  each  locus 
has  one  point  with  the  greatest  slant  range  from  the  launch  point.    This  point  is 
called  the  extreme  optimal  endpoint  (EOF)  for  that  locus.    Since  each  locus  has  an 
EOF,  an  envelope  of  EQP's  is  generated ,  which  represents  the  absolute  maximum 
distance  the  missile  can  reach  before  burnout.    This  envelope  is  shown  as  the  dot- 
dash  cuirve  tangent  to  each  locus  of  optimal  endpoints  in  Figure  14.  .    /     ' 


\  ' 
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Time  History  of  the  Control  Variable 

Figures  7  to  13  present  a  helpful  physical  picture  of  the  optimal  trajectories, 
but  the  desired  quantity  is  the  time  history  of  the  control  variable,  B.    These  histo- 
ries are  presented  in  Figures  15  to  21  ;  they  represent  the  thrust  angle-time  re- 
lationships -which  would  result  in  the  optimum  trajectories  shown  in  Figures  7  to  13. 

Comparison  Trajectories 

In  the  previous  paragraphs,  the  trajectories  obtained  were  referred  to  as 
"optimal"  —  that  is,  minimum-time  paths.  However,  it  must  be  remembered  that 
the  conditions  used  in  deriving  the  trajectory  equations  were  only  necessary  and  not 
sufficient  conditions  for  a  minimum.  Therefore,  it  would  be  instructive  to  compare 
these  trajectories  with  some  arbitrary  trajectories  to  the  same  points.  Due  to  the 
difficulty  of  trying  to  find  arbitrary  trajectories  which  pass  through  given  endpoints, 
this  comparison  is  very  hard  to  make  in  a  numerical  study. 

This  difficulty  can  be  circumvented  in  the  following  way  :    each  locus  of  optimal 
endpoints  represents  a  15-second  constant-time  curve.    If  any  arbitrary  trajectory 
passed  through  the  locus  and  had  an  endpoint  which  fell  outside  the  locus,  then  it 
must  have  crossed  the  locus  at  some  time  earlier  than  15  seconds,  and,  therefore, 
the  "optimal"  trajectory  to  the  point  which  it  crossed  is  not  optimal  at  aU.    Thus,  if 
tile  trajectories  found  in  the  present  study  are  indeed  optimal,  then  all  arbitrary 
trajectories  with  the  same  essential  initial  conditions  should  have  endpoints  which  . 
lie  on  or  within  the  locus  of  optimal  endpoints  corresponding  to  those  initial  condi- 
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tions. 

The  only  comparison  trajectories  used  in  this  study  were  gravity  turn  trajecto- 

ries  (defined  as  trajectories  with  B  S  0) .    For  a  valid  comparison,  the  optimal  and 

comparison  trajectories  must  have  the  same  essential  initial  conditions  —  in  this 

case,  the  same  x,,  h  ,  V   and  A  .    Gravity  turns  for  A^  -  30,  40,  50,  60,  70,  80 
i      i      i  i  i 

and  90  degrees,  respectively,  are  shown  as  the  dot-dash  curves  in  Figures  7  to  13. 
Within  the  limits  of  accuracy  used  in  the  calculations,  the  endpoint  of  each  gravity 
turn  fell  on  the  locus  of  optimal  endpoints  in  each  case.      In  other  words ,  the 
gravity  turns  gave  the  same  performance  as  the  optimal  trajectories. 

In  the  next  section,  an  interception  scheme  is  presented  using  the  concept  of  an 
envelope  of  EOP's  ;  since  gravity  turns  were  found  to  yield  the  same  performance, 
the  same  argument  holds  if  the  locus  of  gravity  turn  endpoints  is  used  in  place  of  the 
envelope  of  EOP's. 

Interception  Scheme 

Once  the  envelope  of  EOP's  has  been  found,  it  can  be  rotated  about  the  h-axls  to 
form  an  "umbrella  of  EOP's"  over  the  launch  site.  The  coordinates  of  each  point  on 
the  umbrella  correspond  to  a  unique  set  of  initial  conditions.  The  time  of  arrival  of 
the  interceptor  at  any  point  if  fired  "now"  is  a  known  quantity.  When  an  enemy  war- 
head is  detected,  its  future  path  can  be  determined,  and  the  point  (in  space  and  time) 
at  which  it  crosses  the  umbrella  can  be  found.    When  the  time  of  crossing  and  the 

For  the  case  of  A    ■  30  degrees, the  gravity  turn  did  not  pass  through  the  locus 
of  optimal  endpoints  ;  tnus  no  direct  comparison  could  be  made. 
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time  of  arrival  Kiatch,  the  interceptor  is  fired  iviih  the  initial  conditions  unique  to 

that  point. 

Conclusions 

It  has  been  shown  that,  within  the  limits  of  accuracy  of  the  calculations,  arbi- 
trary trajectories  yielded  the  same  performance  as  the  optimal  trajectories  in  the 
present  problem.     This  was  thought  to  be  due  to  the  fact  that  the  natural  end  condi- 
tions for  the  variational  problem  could  not  be  satisfied  precisely  and  that  the  tra- 
jectories were,  therefore,  not  true  optimums.    Another  factor  was  the  extremely 
high  thrust  and  correspondingly  high  acceleration,  which  produced  nearly  straigjit- 
line  trajectories,  whether  "optimum"  or  arbitrary.  , 

Using  the  flooding  technique ,  an  "utnbrella"  of  extreme  optimal  endpoints  has 
been  determined  over  the  launch  site ,  with  each  point  on  the  umbrella  correspond- 
ing to  a  unique  set  of  initial  conditions.    An  apparently  feasible  interception  scheme 
using  this  concept  has  been  proposed. 
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APPENDIX  1 


THE  PROBLEMS  OF  BOLZA  AND  MAYER  (5) 


Bolza  Problem 


The  problem  of  Bolza  in  the  Calculus  of  Variations  is  stated  as  follows: 
Given  the  functional, 


){  +       H(t,yk,yk 


I-   (G(t,yj^))t  +     I  H(t,y]j,yij)dtwherek»     1,2,. ..,n, 

find  the  set  of  functions  yi^(t) ,  k  ■  1,  2, . . .  ,n,  which  minimizes  I  while  satisfying 
the  constraints  } 

Pj(t,yj^,yjj)  -0  ,    •       j  -  1,2,...,  p(les8  than  n) 

and  the  boundary  conditions 

Wj.(ti,yj^)  -  0  r  -   l,2,...,q  .  . 

and 

'Wr(tj,yjj^)  -  0  r-  (^l,(jf2, . . .  ,s  (less  than  2n+2), 

Necessary  conditions  for  the  existence  of  an  extremal  for  I  are  that  the  set  of 
yjf's  satisfies  the  Euler- Lagrange  equations,  • 

d*F/d*yjj  -  d/dt(d  F/d  yj^)-    0,  k-    l,2,...,n 

where 

F  »  H  +  IlP. 
5-1  3  3° 

* 
In  the  preceding  two  equations,  d  stands  for  the  symbol  for  partial  differentiation, 

S  stands  for  the  standard  summation  sign  (sigma) ,  and  L^  denotes  a  Lagrange 

J 


la 


Multiplier  which,  in  general,  is  variable  with  time. 

The  set  of  yj^'s  which  extremizes  I  must  also  satisfy  the  p  constraint  equations, 
so  the  differential  system  to  solve  consists  of  n+p  equations  in  the  n  dependent  var- 
iables and  the  p  Lagrange  Multipliers.  ^       . 

The  set  of  yj^'s  must  also  satisfy  the  transversaUty  condition  shown  below  for 
all  dx,  dyjj,  dG  consistent  with  the  boundary  conditions  : 


(dG  +  (F  -  §(d  F/d  yi^)yj^dx+  S(d  F/d  yj^dyj^)  .  -    0, 
k»  1  k»l 


where,  as  before,  d    denotes  a  partial  derivative  and  S  represents  the  summation 
sign. 

Mayer  Problem 
The  Mayer  problem  is  the  special  case  of  the  Bolza  problem  which  occurs 
when  H  "  0  in  the  original  formulation.    For  the  Mayer  problem,  all  of  the  fore- 
going equations  hold.  ' 

In  the  present  problem,  time  is  the  quantity  to  be  minimized,  so  the  function 

f  •'••■.■• 

(G(t,yi^)j    is  merely  tf  -  ti.  - 


■  / 
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APPENDIX  2 


ATMOSPHERIC  DENSITY  PROFILE  AND  DRAG  FUNCTIONC  6) 


The  atmospheric  density  profile  used  in  this  report  was  derived  by  using  the 
basic  differential  equation  for  pressure  vs.  altitude  : 

-dp/p  -  1/R  (dh/T) . 
In  this  equation  p,h,T,  and  R  denote  the  pressure,  altitude,  temperature,  and  gas 
constant.    Plots  of  temperature  vs.  altitude  were  obtained  from  the  NACA  Standard 
Atmosphere,  and  the  above  equation  was  integrated  for  three  altitude  ranges  — 
tropospheric,  stratospheric,  and  above  stratospheric.    The  relation 

Density-    p/gRT 

was  then  used  to  find  the  density  profile.    The  results  were  : 

4.25  ' 

hless  than  35, 000  feet:   Density-     .002378(1  -  .0000069h) 

h  between  35, 000  and  70, 000  feet:    Density-   .  000727  exp  ^^^^°  "  ^ 


h  over  70 ,  000  feet:     Density  > 


21000 

.000136 
((317+.00109h)/  393)^^-^ 


The  units  on  density  and  h  are  slugs  /ft    and  feet,  respectively. 


Drag  Function 


From  Reference  6,  the  drag  equation  is 

D  .  C„(Rho)V^A/2, 
D 

where  Rho,  A,  and  C_.  are  the  atmospheric  density,  frontal  area  of  the  missile, and 


the  drag  coefficient,  respectively.  '  .  ! 

The  profile  drag  coefficient  for  a  streamlined  body  of  revolution  was  given  as 
0.05  ;  since  no  lift  is  present,  this  is  the  total  drag  coefficient.    The  frontal  area       : 
of  the  missile  was  20  square  feet.    Substituting  these  values  into  the  drag  equation,      •  ■ 
we  have  '  '.  - ,  •  . . 

D  .  (Rho)V^/2,     .  .;■  ■  '■:'. 

which  is  the  equation  used  in  this  study. 


.i    .    ■'  ; 
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APPENDIX  3 


» 

RUNGE-KUTTA-GILL  NUMERICAL  INTEGRATION  SCHEME  (4) 


The  Runge-Kutta-Gill  (RKG)  process  is  one  of  many  variations  of  the  basic 

Runge-Kutta  method.    It  is  called  a  "fourth  order"  method  because  the  formulas 

2      3  4 

for  y.    i^,  are  exact  for  fj=  1,  x,  x  ,  x  ,  and  x  . 

Given  a  set  of  differential  equations 

«  ■  .  ■ 

yj  t>  f.(yi'0  »  i  ■  1,2,...  ,N, 
the  domain  of  integration  is  broken  up  into  P  intervals  ;  the  size  of  each  interval 
(step  size)  is  H"    1/P.    If  the  vdlues  of  y.  are  known  at  a  particular  point  in  the 
domain,  say  point  j,  these  values  are  denoted  by  the  notation  y.   ...    The  values  of 
y       .  are  then  given  by  the  following  formulas: 

-1  ^i  ■■-%m :  ■■ '  ■ 

yj^^j^l-yj  .j+l/6(Ai  +  2(l-(l/2)   )Bj+2(H-(l/2)     )Cj  +  Dp  ^ 

where  .    '    ■    '  ■^-     ':  ■•' 

A^  -Hfj(Xj,)^),  W- 

'        Bj    -Hfj(Xj  +  H/2,y  +A/2), 


-1 


.  D,  -  Hf.(x  +  H,  y    -(1/2)     B  +  (1  +  (1/2)     )C) . . 

•      ...         J  J  i  i      .• 


/  - 


7 
/ 


J*   :• 


-1  -1 

Cj   -  Hf  (X  +H/2.y    -  (1/2  -  (1/2)      )A  +(1  -  (1/2)     )B).  ;  '  "l^t     \ 
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APPENDIX  4 


BOUNDED  CONTROL  (^ 


In  this  study,  the  control  variable  was  unbounded.    A  more  realistic  assumption 
would  have  been  that  the  control  variable  (thrust  angle)  has  limits  ---for  example, 
i  10  degrees.    In  order  to  account  for  this,  the  following  constraint  equation  must 
be  imposed  : 

(B+  10)  (10  -  B)  -    K  . 
K  is  a  new  real  variable.    This  is  a  new  holonomic  constraint  on  the  problem,  and 
gives  rise  to  another  Euler- Lagrange  equation  involving  K  and  another  Lagrange 
Multiplier  L5.    This  equation  is  ,      .   . 

KL5-    0.  ..'■■.■',.    (A4. 1) 

The  Euler-Lagrange  equation  for  B  is  also  changed,  becoming 

LgTsinB  -  L4TC0SB  -  L52B  -  0.  (A4.2) 

The  addition  of  these  new  equations  greatly  changes  the  nature  of  the  problem, 
primarily  by  introducing  discontinuous  solutions  for  B(t) .  Due  to  the  limited  scope 
of  this  study,  only  continuous  solutions  with  unbounded  control  were  studied. 
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A  pproved  by: 


The  anti -ballistic  missile  system  which  the  United  States  is  now  deploying  in- 
volves two  missiles,  the  Spartan  and  the  Sprint,  for  high  and  low  altitude  targets 
respectively.    This  report  focuses  on  the  low  altitude  portion  with  the  following  two 
purposes:    to  derive  the  necessary  conditions  for  minimum -time  trajectories  for  a 
low  altitude,  hi^  acceleration  interceptor  such  as  the  Sprint  and  to  integrate  numer- 
ically the  resulting  equations  so  as  to  find  some  optimal  trajectories  for  a  proposed 
model  interceptor  in  the  Sprint  class. 

The  physical  assumptions  that  the  vehicle  motion  was  in  a  vertical  plane,  that 
no  lift  was  present,  and  that  the  flight  path  was  determined  by  the  thrust  angle  alone 
were  made.    Mass  was  assumed  to  be  a  linear  function  of  time,  drag  was  considered 
to  be  a  function  of  altitude  and  velocity,  and  thrust  was  assumed  to  be  constant.    No 
bounds  were  placed  on  the  control  variable  (thrust  angle) . 

The  problem  was  formulated  as  a  Mayer  problem  in  the  Calculus  of  Variations, 
and  the  Euler- Lagrange  equations  and  transversality  condition  were  derived.     The 
number  of  essential  boundary  conditions  allowed  was  determined,  and  natural  end 
conditions  were  found  from  the  transversality  condition. 

A  model  interceptor  was  proposed,  and,  using  data  for  this  model,  the  Euler- 
Lagrange  equations  were  integrated  numerically  to  missile  burnout.    A  locus  of 
"extreme  optimal  endpoints"  was  determined  which  represented  the  locus  of  points 
farthest  from  the  launch  site  at  which  interception  could  occur.    The  concept  was  ex- 
tended to  three  dimensions  by  rotating  this  locus  around  the  altitude  axis,  thus  gen- 
erating an  "umbrella"  of  optimal  endpoints  over  the  launch  site.    Each  point  on  this 


umbrella  corresponded  to  a  unique  set  of  initial  conditions.    When  an  enemy  tra- 
jectory was  found  to  pierce  the  umbrella  at  a  certain  point,  the  correct  firing  angle 
and  time  for  the  interceptor  could  thus  be  found. 

It  was  discovered  that  the  optimum  in  the  present  case  was  very  weak;   arbi- 
trary comparison  trajectories  gave  nearly  equal  performance.    This  was  thought  to 
be  due  to  the  fact  that  the  trajectory  was  extremely  thrust-dominated,  causing  it  to 
approximate  a  straight  line. 
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